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Abstract 

The  temperature  distribution  of  a  thermoelectric  cooler  under  the  influence  of  the  Thomson  effect,  the  Joule  heating, 
the  Fourier’s  heat  conduction,  and  the  radiation  and  convection  heat  transfer  is  derived.  The  influence  of  the  Thomson 
effect  on  the  temperature  profiles,  on  the  fraction  of  the  Joule’s  heat  that  flows  back  to  the  low-temperature  side,  and 
consequently  on  the  maximum  attainable  temperature  difference  and  the  maximum  allowable  heat  load  are  emphasized 
and  explored.  The  results  suggest  that  the  cooling  efficiency  of  a  thermoelectric  cooler  can  be  improved  not  only  by 
increasing  the  figure-of-merit  of  the  thermoelectric  materials  but  also  by  taking  advantage  of  the  Thomson  effect.  A 
possible  development  direction  for  the  thermoelectric  materials  is  thus  given. 
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1.  Introduction 

A  thermoelectric  cooler  pumps  heat,  when  the  elec¬ 
tricity  is  applied,  from  a  low  temperature  reservoir  to 
a  high  temperature  reservoir  through  the  Peltier  effect 
[1].  Compared  to  the  traditional  heat  pumps,  the  ther¬ 
moelectric  cooler  has  the  advantages  of  compact,  high 
reliability,  low  maintenance  fee,  no  vibration  and  easy 
control  (because  of  no  moving  parts),  no  refrigerants 
and  direct  electric  energy  conversion.  The  major  prob¬ 
lem  of  thermoelectric  devices  is  poor  cooling  efficiency. 
This  is  fundamentally  limited  by  the  material  properties 
of  the  n-  and  p-type  semiconductors,  regardless  how  cle- 
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ver  the  device  has  been  engineered.  The  performance  of 
a  thermoelectric  material  is  measured  by  a  dimensionless 
thermoelectric  parameter,  usually  written  as  ZF,  where 
T  is  the  temperature  (usually  the  room  temperature) 
and  Z  characterizes  the  material’s  electric  and  thermal 
transport  properties  and  is  commonly  named  as  the  fig¬ 
ure  of  merit  of  the  material.  The  commonly  known  ther¬ 
moelectric  materials  have  ZT  values  between  about  0.6 
and  1.0.  Different  efforts  have  been  made  to  produce 
semiconductor  materials  with  large  Seebeck  coefficients, 
good  electrical  conductivity,  and  poor  thermal  conduc¬ 
tivity  [2-5]. 

It  is  worthy  to  note  here  that  most  previous  investiga¬ 
tions  [6-8]  of  the  thermoelectric  coolers  considered  only 
the  Fourier’s  heat  conduction  and  the  Joule’s  heating. 
The  radiation  and  convection  heat  transfers  between 
the  semiconductors  and  the  ambient  gas  are  found  not 
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negligible  for  micro-coolers  [7,8].  Little  attention  how¬ 
ever  has  been  paid  to  the  influence  of  the  Thomson  effect 
(charge  carriers  must  release  or  absorb  heat  in  order  to 
overcome  the  non-homogeneity  caused  by  the  tempera¬ 
ture  gradients  in  the  semiconductors)  until  recently 
[9,10].  In  the  present  work,  the  influence  of  the  Thomson 
effect  on  the  performance  of  a  thermoelectric  cooler, 
such  as  its  maximum  attainable  temperature  difference 
and  maximum  allowable  heat  load,  is  of  interest  and 
investigated.  These  results  provide  the  theoretical  bases 
for  the  possible  development  of  thermoelectric  materials 
and  for  an  optimal  design  of  a  thermoelectric  cooler. 


2.  Thermal  model 

The  basic  unit  of  a  thermoelectric  cooler  (TE  cooler) 
is  a  thermocouple  composed  primarily  of  an  n-type  and 
a  p-type  semiconductor  element  placed  electrically  in 
series  and  thermally  in  parallel,  as  shown  in  Fig.  1. 
The  attempt  is  to  pump  heat  from  the  cold  side  (x  =  0) 
through  the  Peltier  heat,  transfer  it  to  the  hot  side 
(x  =  L )  by  current,  and  finally  release  it  to  the  ambient 
by  some  external  cooling  device.  The  cooling  power 
however  cannot  be  as  large  as  the  Peltier  heat  pumping 
rate,  because  additional  heat  sources  such  as  the  Joule’s 
dissipation  heat  and  the  Thomson  heat  are  generated 
and  heat  will  diffuse  back  to  the  cold  side  through  the 
Fourier’s  heat  conduction.  In  the  present  study,  a  one¬ 
dimensional  analysis  will  be  employed,  which  is  proper 
when  the  length  of  the  beam  is  sufficiently  larger  than 
its  width  and  thickness.  According  to  the  nonequilib¬ 
rium  thermodynamics,  the  steady  temperature  distribu¬ 
tion  T(x)  is  governed  by 

o-Uudi)+‘¥AT^‘j-TA',-yp(T-T") 

(i) 

Terms  on  the  right-hand-side  represent  respectively 
the  net  axial  thermal  conduction  heat  transfer  rate,  the 
Joule’s  heat  generation  rate,  the  Thomson  heat  genera¬ 
tion  rate,  and  the  thermal  radiation  and  convection  heat 


transfer  rates.  Symbols  2  and  a  represent  the  thermal 
and  electrical  conductivities;  ee  and  je  are  the  electron 
charge  and  flux  (therefore  the  current  is  I  =  Aee  je);  (3 
is  the  Thomson  coefficient;  A  and  P  are  the  local 
cross-sectional  area  and  perimeter  of  the  semiconductor 
element.  The  sign  of  the  term  associated  with  the  Thom¬ 
son  heat  depends  on  the  current  direction.  It  is  a  minus  if 
the  current  flows  in  the  positive  x-direction  (the  p-type 
element  in  Fig.  1);  and  a  plus  otherwise  (n-type).  Finally, 
the  radiation  heat  transfer  rate  has  been  linearized  under 
the  assumption  of  small  temperature  difference  and  com¬ 
bined  with  the  convection  heat  transfer  coefficient  h  such 
that  y  =  4saB T3^  +  /z,  where  g,  crB,  and  are  the  emis- 
sivity,  the  Stefan-Boltzman  constant,  and  the  ambient 
gas  temperature. 

Further  simplifications  can  be  made  by  assuming 
constant  properties  (A,  (3,  and  a),  constant  A  and  P, 
and  constant  current  /.  Emphasized  is  the  assumption 
of  constant  Thomson  coefficient.  By  definition,  the 
Thomson  coefficient  /3  is  related  to  the  Seebeck  coeffi¬ 
cient  oc  by 


(2) 


where  the  subscription  “p”  represents  a  constant  pres¬ 
sure  condition  under  which  the  derivative  is  taken. 
Therefore,  a  constant  Seebeck  coefficient  means  an  ab¬ 
sence  of  the  Thomson  effect.  On  the  other  hand,  a  con¬ 
stant  Thomson  coefficient  implies  a  Seebeck  coefficient 
having  a  logarithm  dependence  on  the  temperature. 


2.1.  Thomson  effect  on  the  temperature  profiles 


With  the  boundary  conditions  T(0)  =  Tc  and 
T(L)  =  Th,  the  solution  of  Eq.  (1)  is 
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Fig.  1.  Schematic  diagram  of  a  thermoelectric  cooler. 
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where  x*  =  x/L,  ©  =  (T-T^/AT,  ©c  =  (T^T^/AT, 
A T=  Th-Tc,  P*  =  PyL/K,  and  Y  =  I2R/KAT.  Symbols 
R  =  LI  o  A  and  K=  2A/L  are  the  electric  resistance  and 
thermal  conductance  of  thermoelement.  The  dimension¬ 
less  parameters  Ki  and  k2  are  the  two  distinct  real  roots 
of 


0  =  k2  =f  £k  -  P* 


(4) 


with  £  =  fllK  representing  the  ratio  of  Thomson  heat  to 
the  conduction  heat.  Noticed  is  that  the  parameters  K\ 
and  k2  are  now  both  functions  of  the  operating  current 
/  because  of  the  Thomson  effect. 

To  highlight  the  influence  of  the  Thomson  effect  on 
the  temperature  distribution,  we  choose  y  =  0  for  the 
time  being.  The  two  roots  of  Eq.  (4)  are  then  ±£  and  0 
and  the  temperature  distribution  becomes 


T(x)-Te  1  ~  exp  (±Qf)  x 

Tb-Tc  1  +U'  l-exp(±5)  T 


(5) 


for  the  p-type/n-type  element  with  ^=I2  Rip  I  AT  being  the 
ratio  of  Joule’s  heat  to  the  Thomson  heat.  Noticed  is 
that  two  influence  dimensionless  parameters,  £  and  f, 
appear  separately,  unlike  the  case  when  the  Thomson 
effect  is  ignored  and  the  temperature  profile  is  then 


n*)-rc  =  i  Mx_J\  x 

Th-Tc  2a\L  L2J+L 


(6) 


where  £ f  is  the  ratio  of  the  Joule’s  heat  to  the  con¬ 

duction  heat.  Shown  in  Fig.  2  are  temperature  profiles  of 
the  p-type  element  under  several  values  of  the  Thomson 
coefficients.  As  seen,  when  the  Thomson  coefficient  is 
negative  (pp  <  0),  the  temperature  near  the  cold  side  is 
increased  by  the  Thomson  heat  released  by  the  charge 
carriers.  This  temperature  rise  will  make  a  reduction  in 
the  cooling  power  of  the  TE  cooler.  On  the  other  hand, 
if  Pp  >  0,  the  temperature  gradient  is  everywhere  positive 
and  charge  carriers  absorb  heat  all  the  way  as  the  cur- 


1.2 
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Fig.  2.  The  influence  of  Thomson  effect  on  the  temperature 
profiles  where  T  =  2  and  -6  ^  ^  6  are  chosen. 


rent  moves  from  x  =  0  to  x-L.  This  implies  that  most 
of  the  Joule’s  heat  and  the  Fourier’s  conduction  heat 
are  luckily  absorbed  and  carried  to  the  hot  side  by  the 
charge  carriers.  Therefore,  for  the  sake  of  improving 
the  cooling  power,  a  TE  cooler  should  employ  a  p-type 
(n-type)  semiconductor  material  with  a  positively  (nega¬ 
tively)  large  Thomson  coefficient,  in  addition  to  a  low 
electric  resistance  and  a  low  thermal  conductance. 

2.2.  Thomson  effect  on  the  cooling  power 

The  cooling  power  ( qc )  is  the  net  heat  flow  rate  out  of 
the  cold  side  and  is  computed  as 

<7C  =  «pn {Te)rre  -  yp  ^  (0)  -  Vfn  ^  (0)  (7) 

where  apN  =  ap  —  an,  and  subscriptions  p  and  n  repre¬ 
sent  p-type  and  n-type  materials  respectively.  Using 
Eq.  (3),  it  can  be  shown  that  the  cooling  power  is  equal 
to 

qe  =  «pN(rc)/rc  -  K(Th  -  Tc)  -  H(TX  -  Te)  -  I2R, 


(8) 

where 

K  =  Kpf  {kp\  ,  Kp2 )  T  Kn f  (xn i ,  Kn 2) 

(9) 

H  =  Kpg(KpX ,  JCp2)  +  Kng(KnX ,  Knl) 

(10) 

R  =  Rph{Kp\:  Kp2 )  T  Rnh(Kn\:  Kn 2) 

(11) 

and  functions  /,  g,  and  h  are  defined  as 

/(/Cl,  k2)  =  — - — 

J  v  ’  eKl  —  eK2 

(12) 

Kie*2  -  K2eKl 

g(KUK2)=  eK1_gK2  f(KUK2) 

(13) 

K\K2 

(14) 

Firstly  of  all,  it  is  noticed  that  the  fraction  of  Joule’s 
heat  flowing  to  the  cold  side  is  no  longer  one  half  as 
usual  when  no  Thomson  effect  is  concerned.  Instead,  it 
is  determined  by  the  function  h(Ki,K2)  which  ranges  be¬ 
tween  0  and  1.  The  influence  of  the  radiation  and  con¬ 
vection  heat  transfer  lies  not  only  in  H  but  also  in  K 
and  R  because  K\  and  k2  are  also  functions  of  y  as  ob¬ 
served  from  Eq.  (4). 

Again,  to  illuminate  the  Thomson  effect,  we  stick  to 
cases  in  which  y  =  0.  In  those  cases,  H  =  0  and  the 
so-called  herein  modified  thermal  conductance  K  and 
the  modified  electric  resistance  R  become 


K=KptlK(Qp)+Kn>lK(-Q 


(15) 


R  ^p7/?(Cp)  “H  On) 


(16) 
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Fig.  3.  The  dependence  of  rjK  and  rjR  on  the  dimensionless 
operating  current 


where 

[~1_f  as<^<l  (17) 

=T^  +  \R1\~  J2  as^<L  (18) 

A  value  of  rjK  greater  than  one  implies  the  additional 
(Thomson)  heat  that  must  diffuse  to  the  cold  side 
through  the  Fourier’s  conduction.  The  value  of  rjR 
(0  ^  Or  ^  1)  on  the  other  hand  tells  the  fraction  of 
Joule’s  heat  that  will  flow  to  the  cold  side.  Fig.  3  shows 
the  dependence  of  rjK  and  ^  on  the  operating  current 
nondimensionalized  in  terms  of  the  thermal  conductance 
and  the  Thomson  coefficient.  It  is  seen  that  both  the 
modified  thermal  conductance  and  the  modified  electric 
resistance  can  be  significantly  reduced  by  properly 
choosing  semiconductor  materials  with  f3v  >  0  and 

Pn<0. 


3.  Performance 

Given  a  temperature  difference  A  T-  Th  —  Tc,  there 
may  exist  some  optimum  operating  current  that  maxi¬ 
mizes  qc.  The  maximum  attainable  temperature  differ¬ 
ence  A  Tmax  is  the  temperature  difference  when  the 
corresponding  maximum  qc  is  zero.  And  the  maximum 
allowable  heat  load  Amax  is  defined  as  the  maximum 
cooling  power  when  AT  =  0.  The  involved  algebra  is  a 
little  troublesome  because  the  modified  thermal  con¬ 
ductance  and  modified  electric  resistance  are  now  func¬ 
tions  of  I.  For  simplicity  and  for  the  sake  of 
illuminating  the  Thomson  effect,  y  =  0  is  assumed  again. 
Also  assumed  is  that  the  p-type  and  n-type  materials 


have  the  same  constant  thermoelectric  properties  except 
that  apn(rc)  =  2a(71c)  >  0  and  Pp  =  -/?n  =  P  >  0. 

3.1.  Maximum  attainable  temperature  difference 

After  some  algebra,  it  can  be  shown  that  the  maxi¬ 
mum  attainable  temperature  difference  A  Tmax  and  the 
corresponding  optimum  operating  current  /opt.T  are 
determined  by 

ZTC  =  ^(1  —  e-{optT)  (19) 

ZATmax  =  j  (£opt.T  +  e-«~T  -  1 )  (20) 

where  £opt.T  =  piopt.T!K  and  Z  =  gol2IX  is  the  Figure-of- 
Merit  evaluated  at  Tc.  The  influence  of  the  Thomson 
effect  can  be  explained  in  the  following  way.  Given  Tc 
and  a  ( Tc ),  Eq.  (19)  gives  the  /^-dependent  optimum 
dimensionless  operating  current  and  Eq.  (20)  in  turn 
gives  the  maximum  attainable  temperature  difference. 
Moreover,  from  Eq.  (19),  it  is  seen  that  when  ZTC  is  less 
than  the  value  of  dp  (or  Tc  is  less  than  XI  g a/?),  there  exists 
an  £opt.T  with  a  corresponding  maximum  attain¬ 
able  temperature  difference.  Otherwise,  there  exists  no 
so-called  optimum  operating  current  and  there  is  no  limit 
on  the  temperature  difference.  This  is  because  most  of  the 
Joule’s  heat  will  just  be  absorbed  as  Thomson  heat  and 
carried  to  the  hot  side  by  the  current,  becoming  not  an 
obstacle  in  creating  the  desired  temperature  difference. 
All  of  the  above  are  illustrated  in  Fig.  4.  In  remark,  both 
parameters,  ZTC  and  ct/fi,  play  roles  in  determining  the 
performance  of  a  thermoelectric  cooler. 


ZTC 


Fig.  4.  The  maximum  attainable  temperature  difference 
ZATmax  versus  ZTC.  Values  of  dp  range  in  a  geometric 
progression  with  a  ratio  of  1.6  starting  from  0.2.  The  dotted 
curve  is  the  one  corresponding  to  P  =  0  (or  dp  =  oo). 
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3.2.  Maximum  allowable  heat  load 


Using  Eqs.  (8)  and  (15)-(18),  the  maximum  allowable 
heat  load  Amax  can  be  found  to  be 


a2  max 


2£opt.N  L  ^opt.Ne^  \ 

( 1  —  q^°pln  )  y  ( 1  —  )  J 


(21) 


P  rr  _  1  ,  2l=opt.N  , 

—  h  -  1  \  - e - 1 - S 

a  1  —  ec°ptN  (1  —  e^pt-N) 


(22) 


where  N*mllx  =  ZNmm/K  and  copt,N  =  P  Iovt.vslK  (the  opti- 
mum  dimensionless  operating  current  that  results  in  a 
maximum  allowable  heat  load),  and  is  shown  in  Fig.  5. 
A  careful  examination  of  Eq.  (21)  finds  that  there  exist 
two  critical  values  of  ZTh/3l a,  namely  ZTh/3l a  =  1  and 
ZTh/3/oc  ^  1.18.  When  ZTh/3/ot  <  1,  there  exists  an  opti¬ 
mum  operating  current  £opt.N  and  a  corresponding  max¬ 
imum  heat  load  A^ax.  As  1  <  ZTh/3/ a  <  1.18,  there  is  a 
local  maximum  heat  load  at  some  operating  current 
and  a  local  minimum  heat  load  at  some  larger  operating 
current.  As  the  operating  current  becomes  even  larger, 
the  allowable  heat  load  increases  monotonically.  These 
local  extremes  are  presented  by  the  two  branches 
in  the  region  of  1  <  ZThf/oi  <  1.18  in  Fig.  5.  When 
Z7y?/a>1.18,  the  heat  load  increases  monotonically 
with  the  operating  current  and  there  exists  thus  no  so- 
called  maximum  allowable  heat  load.  All  of  the  above 
situations  are  shown  in  Fig.  6,  in  which  the  dimension¬ 
less  heat  load,  N*  =  ZN/K ,  is  plotted  as  a  function  of 
the  dimensionless  operating  current  under  the  constraint 
AT  =  0  for  several  values  of  ZThfl a. 

It  is  not  so  obvious  however  why  the  Thomson  effect 
can  enhance  the  cooling  power  since  the  total  Thomson 
heat  carried  away  by  the  current  is  in  fact  zero 
(PIAT  =  0).  Recalled  is  that  the  Thomson  effect  can 


5 


Fig.  6.  The  allowable  heat  load  under  the  constraint  AT  =  0  as 
a  function  of  the  operating  current  for  various  values  of 
ZTh  p/oL. 


change  the  temperature  profiles  (as  investigated  in  Sec¬ 
tion  2.1)  and  therefore  direct  a  major  portion  of  the 
Joule’s  heat  flowing  to  the  hot  side  (in  the  sense  that 
charge  carriers  absorb  heat  near  the  cold  side  and  release 
the  same  amount  of  heat  in  the  neighborhood  of  the  hot 
side).  The  cooling  power  is  thus  increased.  When  the 
Thomson  effect  is  sufficiently  strong  ( ZThf/oi  >1),  the 
fraction  of  Joule’s  heat  and  the  Fourier’s  conduction 
heat  flowing  to  the  cold  side  become  negligible  and  con¬ 
sequently  the  cooling  power  (or  the  allowable  heat  load) 
can  be  made  almost  equal  to  the  Peltier  heat,  which  is  in 
turn  proportional  to  the  operating  current.  This  limiting 
behavior  however  can  be  only  slowly  approached  as  seen 
from  Fig.  3  or  from  Eq.  (18)  (i.e.  rjR  «  l/£  as  £  >  1.) 


Fig.  5.  The  extreme  values  of  allowable  heat  load  under  the 
constraint  AT  =  0  versus  ZTh/3/ot.  Dotted  line  is  the  result 
without  considering  Thomson  effect. 


4.  Illustration 

At  the  end  of  this  paper,  we  illustrate  the  present 
analysis  by  an  example.  The  properties  of  the  p-type 
and  n-type  elements  are  so  chosen  [7,8]  that,  op  = 
<rn  =  0.1  (pQm)-1,  ap  =  2n  =  1.6W/mK,  ap  =  —  an  =  a, 
and  finally  fp  =  —fn  =  f  >  0.  The  geometry  is  chosen 
to  be  Lp  -  Ln  -  1  mm,  Pp  =  Pn  =  0.4mm,  and  Ap  =  An  = 
0.01mm2  The  radiation-convection  heat  transfer  coeffi¬ 
cient  is  y  =  50W/m2K  (a  reasonable  value  for  air  at  nor¬ 
mal  atmosphere).  The  ambient  temperature  T^  is  taken 
to  be  300  K. 

Fig.  7  shows  the  temperature  difference  as  a  function 
of  the  operating  current  when  no  heat  load  is  applied. 
The  cold-side  temperature  is  fixed  at  250  K  and 
a (Tc)  =  185  pV/K  is  assumed.  As  seen,  when  the  Thom¬ 
son  effect  is  ignored,  the  optimum  current  is  about 
0.051  A  and  the  maximum  temperature  difference  is 
about  55  K.  The  value  of  2ATmax/rc  is  about  0.44,  which 
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Fig.  7.  The  attained  temperature  difference  versus  the  operating 
current  under  no  heat  load.  The  Thomson  coefficients  are  given 
in  pV/K. 


the  Seebeck  coefficient  of  a  one-band  (parabolic  band) 
bulk  material  is  determined  by 


( $F_ 3/2 


(23) 


where  f*  =  f*(7)  is  the  Fermi  energy  nondimensionalized 
by  kBT  and 


xldx 

exp(x  -  C)  +  1 


(24) 


Equation  (23)  suggests  that  it  may  be  possible  to  con¬ 
trol  the  temperature  dependence  (derivative)  of  the  See¬ 
beck  coefficient  (and  thus  control  the  Thomson 
coefficient)  through  a  control  of  the  Fermi  energy  and 
even  its  temperature  dependence,  say  by  doping. 
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is  smaller  than  ZTC& 0.53  because  of  the  convection  and 
radiation  effect.  The  maximum  attainable  temperature 
difference  is  seen  to  increase  with  increasing  Thomson 
coefficient.  It  is  about  114K  (2ATmax/rc^0.91)  when 
/?  =  200pV/K.  Finally,  the  temperature  difference  be¬ 
comes  unlimited  when  /3  is  slightly  larger  than  the  critical 
value  of  2looiTc  «  346pV/K. 


5.  Conclusion 

A  one-dimensional  thermal  analysis  for  the  perform¬ 
ance  of  thermoelectric  cooler  has  been  conducted  under 
the  influence  of  the  Thomson  effect,  the  Joule  heating, 
the  Fourier’s  heat  conduction,  and  the  radiation  and 
convection  heat  transfer.  The  temperature  distribution 
and  the  cooling  power  are  analytically  derived.  For 
those  cases  that  the  Thomson  effect  is  important,  two 
dimensionless  influence  parameters,  i.e.,  the  ratio  of 
the  Thomson  heat  to  the  conduction  heat  and  the  ratio 
of  Joule’s  heat  to  the  Thomson  heat,  play  roles.  Both  of 
the  temperature  distribution  and  the  cooling  power  are 
affected  by  the  Thomson  effect.  The  fraction  of  both 
Fourier’s  conduction  heat  and  Joule’s  heat  that  flows 
back  to  the  cold  side  can  be  significantly  reduced.  It  re¬ 
sults  in  higher  values  of  maximum  attainable  tempera¬ 
ture  difference  and  the  maximum  allowable  heat  load, 
providing  that  the  thermoelectric  cooler  has  p-type  (n- 
type)  semiconductor  material  with  the  positive  (nega¬ 
tive)  Thomson  coefficient. 

Although  a  material  with  a  significant  Thomson  ef¬ 
fect  is  still  limited  due  to  the  material  capability  at  the 
present  time,  however  this  paper  has  highlighted  a  pos¬ 
sible  development  direction  for  the  thermoelectric  mate¬ 
rials  in  the  future.  For  example,  it  was  shown  [11]  that 
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